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Abstract 

The study of the transition between galactic and extragalactic cosmic rays can shed more light on the end of the Galactic cosmic rays 
spectrum and the beginning of the extragalactic one. Three models of transition are discussed: ankle, dip and mixed composition 
models. All these models describe the transition as an intersection of a steep galactic component with a flat extragalactic one. 
Severe bounds on these models are provided by the Standard Model of Galactic Cosmic Rays according to which the maximum 
acceleration energy for Iron nuclei is of the order of E^ x ~ 1 x 10 17 eV. In the ankle model the transition is assumed at the 
ankle, a flat feature in the all particle spectrum which observationally starts at energy E a ~ (3 — 4) x 10 18 eV. This model needs 
a new high energy galactic component with maximum energy about two orders of magnitude above that of the Standard Model. 
The origin of such component is discussed. As observations are concerned there are two signatures of the transition: change of 
energy spectra and mass composition. In all models a heavy galactic component is changed at the transition to a lighter or proton 
component. As a result the ankle model predicts a galactic Iron component at E < 5 x 10 18 eV, while both HiRes and Auger 
data show that at (2 — 5) x 10 18 eV primaries are protons, or at least light nuclei. In the dip model the transition occurs at the 
second knee observed at energy (4 — 7) x 10 17 eV and is characterized by a sharp change of mass composition from galactic Iron 
to extragalactic protons. The ankle in this model appears automatically as a part of the e + e~ pair-production dip. The mixed 
composition model describes transition at E ~ 3 x 10 18 eV with mass composition changing from the galactic Iron to extragalactic 
mixed composition of different nuclei. In most mixed composition models the spectrum is proton-dominated and it better fits 
HiRes than Auger data. The latter show a steadily heavier mass composition with increasing energy, and we discuss the models 
which explain it. 

Keywords: ultrahigh energy cosmic rays, galactic cosmic rays, cosmic ray theory, cosmic ray experiment 
PACS: 95.85.Ry, 96.40.-z, 95.85.Ry, 98.70.Sa 



X 



1. Introduction 

A clear understanding of the provenance of Cosmic Rays 
(CR) at different energetic regimes, namely their galac- 
tic or extragalactic origin, is of paramount importance in 
unveiling the nature of the sources of these particles, in 
particular at the highest energies. 

The observed energy spectrum of CR has an approxi- 
mately power- law behavior for 11 orders of magnitude in 
energy with several features that can be linked with parti- 
cles propagation and acceleration. This power-law behav- 
ior is most probably indicative of a power-law acceleration 
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spectra, while spectral features may be assigned to changes 
in the origin of particles, their propagation and accelera- 
tion. 

We shall consider these features starting from the high- 
est energies. 

If primary particles are protons, their spectrum must 
show a steepening which begins at Eqzk ~ 50 EcV. It 
is the famous Greisen [l|, Zatsepin and Kuzmin ^ cut- 
off, referred to as GZK, caused by pion photo- production 
energy loss in collisions of protons with the Cosmic Mi- 
crowave Background (CMB) photons. Most probably this 
cutoff is already observed in the HiRes (High Resolution 
Fly's Eye) [H and TA (Telescope Array) Q data, while the 
spectrum steepening observed by the PAO (Pierre Auger 
Observatory) [5[ does not agree well, according to our cal- 
culations, with the predicted GZK shape and position. 

An alternative to the protons as primaries is given by 
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nuclei. Propagating through astrophysical backgrounds 
nuclei lose their energy in photo-disintegration and pair 
production processes, and also due to adiabatic expansion 
of the universe. The steepening of nuclei spectra due to 
photo-disintegration in the interactions with the CMB has 
been shortly mentioned by Greisen [l[ and later consid- 
ered quantitatively in @, 0, B E3| • The beginning of the 
steepening in this case (nuclei in the CMB radiation only) 
corresponds to the intersection of e~e + and adiabatic en- 
ergy losses curves 0, HI ■ The beginning of this steepening 
is lower than in the case of GZK cutoff and the shape of 
this feature differs from that of GZK. Recent accurate cal- 
culations show that the PAO spectrum steepening may be 
explained not only by nuclei interactions. It may be also 
produced as a result of a decreasing acceleration efficiency 
in the vicinity of the maximum energy that sources can 
provide (E max ) (for reviews see The HiRes and TA 



data, in favor of proton composition, and the PAO data, 
in favor of nuclei composition, are further discussed in sec- 
tions 13.21 and 13.31 

There is another feature in the spectrum of UHE ex- 
tragalactic protons, the pair-production dip [l2j |. which, as 
in case of GZK cutoff, can be directly linked to the in- 
teraction of CRs with the CMB. This dip arises due to 
electron-positron production energy loss by extragalactic 
protons interacting with the CMB photons: 



V + 7CMB 



p 



(1) 



This feature has been studied in [13j, [lj, [15|, [16[ . The 
dip was observed with a very good statistical significance, 
X 2 /d.o.f. - 1, by the Fly's Eye, Yakutsk, Akcno-AGASA 
and HiRes detectors , but it is absent in the Auger data 
(see section l3~3l for a more detailed analysis). 

The pair-production dip and GZK cutoff are signatures 
of protons. A confirmation of the shape of these features 
can be considered as an indirect evidence for a proton- 
dominated composition of primary CRs. For nuclei as 
primaries the shapes of the dip and cutoff are strongly 
different. 

A different explanation of the dip has been proposed by 
Hill and Schramm fl8j . They interpreted the dip observed 
in 1980s in terms of a two-component model; the low en- 
ergy component was either galactic or produced by Local 
Supcrclustcr. A similar model was later considered in (l9j . 
The Hill-Schramm's interpretation is widely adopted now. 

From 1960s there was observed a flattening in the Ultra 
High Energy Cosmic Rays (UHECR) spectrum, which was 
called the ankle. Discovery of this feature at the Volcano 
Ranch detector [20| was interpreted as the transition from 
a steep galactic component to a flatter extragalactic one. 
The ankle is clearly seen in the PAO data at E° hB w 4 EcV 



211 ] . In [22| - [23J the ankle is considered as a feature where 
the transition from the galactic to extragalactic cosmic 
rays occurs. In the dip model the ankle appears 

automatically as a part of the pair-production dip. 

Below the ankle, at -E s kn ~ (0.4 — 0.7) EeV, there is a 
faint feature in the spectrum [24| called the second knee. 



It is seen in many experiments (for a review see [2q ) . This 
feature is also often interpreted as the place of transition 
from galactic to extragalactic CRs, with the dip model 
14l . Il5l . Il6| being among these works. 



One more well known feature, the proton knee, occurs at 
E p ~ 4 PcV. It was discovered by the MSU (Moscow State 
University) array in 1958 j2f|. The knee is most probably 
explained by the maximum acceleration energy for protons 
in the Galaxy. In the framework of the rigidity-dependent 
acceleration model it was predicted that the Iron knee is 
located at Zp c = 26 times higher energy, Ep c ~ 0.1 EeV. 
This feature was confirmed by the KASCADE Grande (27l |. 
so that the end of the Galactic CR spectrum is naturally 
connected with the Iron knee. 

Based on the observational features described above, 
three models of the transition from galactic to extragalac- 
tic cosmic rays have been proposed. 

In the ankle models [22| - [23| it is assumed that the tran- 
sition occurs at the flat part of the observed spectrum in 
the energy interval £'* ran8 ~ (3 — 10) EcV. The transition 
energy is given by the intersection of a flat extragalactic 
spectrum and a very steep galactic one. In the majority of 
ankle models the extragalactic component is assumed to 
be pure proton, while the galactic one should be naturally 
represented by Iron nuclei at energies above the Iron knee. 
These models predict a transition from an Iron-dominated 
composition to a proton-dominated one at the ankle en- 
ergy. 

In the mixed composition model [28| it is assumed that 
the extragalactic component consists of nuclei of various 
types. Thus transition here occurs from Iron to lighter 
nuclei of mixed composition; it can occur at the ankle or 
nearby it. 

In the dip model the transition begins at the second 
knee and is completed at the beginning of the dip, at E ps 
1 EeV. The ankle in this model appears as an intrinsic 
part of the dip. Like in the ankle model, the transition 
here also occurs as an intersection of the flat extragalactic 
component (this flatness is especially prominent in the case 
of diffusive propagation) with a steep Galactic spectrum. 
In contrast to the ankle and mixed composition models, 
the dip model predicts an almost pure proton composition 
above E k, I EcV and a pure Iron composition below this 
energy. 

In this review we shall discuss recent modifications of 
these three basic models. The paper is organized as fol- 
lows: in section[5]wc briefly discuss the standard model of 
Galactic cosmic rays, focusing in particular on the end of 
the spectrum, in section|2]wc discuss the spectra of cosmic 
rays at ultra high energy, focusing on their extragalactic 
origin, in section [4] we analyze the models of the transi- 
tion from Galactic to extragalactic cosmic rays. Finally, 
in section [5] we shall present our conclusions. 
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2. The end of the Galactic CRs 

For a long time it was assumed that Supernova (SN) 
explosions may provide the observed cosmic rays. Baade 
and Zwicky [29[ in 1934 were the first who understood that 
SNe energy release can produce enough energy to fill the 
Milky Way with the observed cosmic rays. Arguing with 
the adherents of extragalactic origin of the observed cos- 
mic rays, Ginzburg and Syrovatsky put forward a model 
[30I ] of Galactic Cosmic Rays (GCR). Its basic elements 
were particle acceleration by SN shocks followed by a dif- 
fusive propagation in the Galactic magnetic fields, until 
they leave the Galaxy. 

2.1. Standard Model for Galactic Cosmic Rays (SM GCR) 

In spite of some small disturbing contradictions that will 
be discussed later, one may claim now that the Standard 
Model for Galactic Cosmic Rays exists (see reviews 
32, Ull). It is based on the Supernova Remnant (SNR) 
paradigm and includes four basic elements: 

1. supernova remnants (SNRs) as sources; 

2. SNR shock acceleration; 

3. rigidity-dependent injection as a mechanism providing 
the observed CR mass composition, and 

4. diffusive propagation of CRs in the Galactic magnetic 
fields. 

Let us discuss these elements in more detail. 

1. SNRs are able to provide the observed CR energy 
production in the Galaxy, which can be found [3l| as 



Q cr « CJ CT cM g /x CT , 



(2) 



where w cr ps 0.5 eV/cm 3 is the observed CR energy density, 
c is the velocity of CR particle, M g « 5 x 10 42 g is the total 
mass of galactic gas, and x cr « 7 g/cm 2 is the grammagc 
traversed by CR before escaping the Galaxy. Using these 
numbers one obtains Q« 2 x 10 40 erg/s, which is less than 
10% of the energy release in the form of kinetic energy of 
the SN ejecta per unit time. 

2. The most efficient mechanism of diffusive shock accel- 
eration (DSA) in SNRs was discovered in 1977 - 1978 Q. 
A great progress has been reached during the last decade 
in the theory of acceleration. It was shown that cosmic ray 
streaming instability strongly amplifies magnetic fields up- 
stream. It creates a highlyturbulent field with strength 



up to SB r~i B r~i 10 4 G [35( (for recent works see |36|). 



At each moment of the shock propagation only particles 
accelerated up to maximum energy -Eniax can escape the 
acceleration region {37}. As a byproduct, this solves the 
problem of adiabatic energy losses, since particles escaping 
with E max do not dwell long in the expanding SN shell. 

fi max reaches the highest value at the beginning of the 
Scdov phase and then diminishes due to the shock de- 
celeration. At each moment t the spectrum of escaping 
particles has a narrow peak at E max (t), but the spectrum 



integrated over time acquires a classical E~ 2 shape with a 
flattening at highest energies [371 ] . The maximum acceler- 
ation energy estimated in the Bohm regime of diffusion is 
given by 



B 



W51 



10" 4 G \n g /cm 



2/5 



PcV, 



(3) 



where B is the amplified magnetic field, W51 is the kinetic 
energy of the shell in units 10 51 erg, n g is the upstream 
density of the gas and Z is the charge number of acceler- 
ated nuclei. Thus for protons and Iron nuclei maximum 
energies are 



£™ = 4xS 

771 max 

E 



4 PeV, 
0.1 x B_ 4 EeV. 



(4) 



Here £"™ ax describes well the position of the proton knee 
and E™^* predicts the position of the Iron one. 

An important result on the structure of the knee was 
obtained in the KASCADE-Grande experiment [23 • The 
spectra were measured separately for electron-rich show- 
ers, initiated by protons and light nuclei, and electron- 
poor showers, corresponding to the heaviest nuclei of the 
Iron group. The steepening in the spectrum of electron- 
poor showers was found at E ss 80 PeV. With electron- 
rich showers knee at (3 — 5) PeV, the heaviest nuclei knee 
corresponds to Iron nuclei, as expected in the rigidity- 
dependent acceleration model. 

3. Observations show that nuclei in cosmic rays are sys- 
tematically more abundant than in the interstellar medium 
of the solar neighborhood [H, [3^, 140] : this occurs due 
to the injection of particles in the regime of accelera- 
tion [H, |H, It can be illustrated by the following 
simple consideration flit ]. A particle i from downstream 
(t = A, p) can cross the shock and thus may be injected in 
the regime of acceleration if its Larmor radius Tl{jp) > d, 
where d is the thickness of the shock front. It readily gives 
the relation between a nucleus and proton injection mo- 
menta 



V 



ZeBd/c = Zp? n . 



(5) 



rate for nuclei. 



Therefore, < vf n ^, which provides a higher injection 



The same result may be also obtained in a more formal 
way. Consider a flux of accelerated particles i, Ji{p) = 
Ki(p/p\ n ^)~ lg . Normalizing Ji(p) by the condition 



-in 
c 



Ji(p)dp = r\in u 



(6) 



where n 2 ; is the density of gas i and r\i is the fraction of 
this density injected into acceleration process, one obtains 
a ratio of nuclei and proton CR fluxes 



Ja(p) 
Mp) 



Z 1 ' 



TjpTlp 



(7) 
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Thus, the fraction of nuclei is enhanced by a factor 
Zva/Vp- For numerical calculations of CR nuclei abun- 
dances see [11 III, [II. 



4. CRs propagate in the Galaxy diffusively; they scat- 
ter off small-scale magnetic turbulence which may be de- 
scribed as a superposition of MHD waves with different 
amplitudes and random phases. This process is considered 
(see e.g. 3l|) in the resonance approximation, when the 
particle giro-frequency is equal to the wave frequency in 
the system at rest with the particle motion along the aver- 
age magnetic field. The magnetic field may be considered 
as a combination of an average (constant) field, Bq, and 
a fluctuating component, B. In [3l| the parallel diffusion 
coefficient D\\(E) is calculated assuming D±(E) <C D\\{E) 
(see however numerical simulations [41( which do not sup- 
port this assumption at the highest energies). 

Proton and nuclei spectra, as calculated by Berczhko 
and Volk [H] within SM GCR, are shown in Fig. Q] in 
comparison with data. An agreement with observations 
is quite good at low energies, and the knee is confirmed 
at -Ekn ~ 3 PcV in proton and all-particle spectra. The 
Iron knee located at Ep e ps 80 PeV is the most impor- 
tant prediction of the Standard Model; this feature was 
confirmed by the KASCADE-Grande detector ^ . The 
spectra beyond knees are expected to be very steep. 

2.2. Diffusive propagation: qualitative estimates 

Diffusive propagation is a phenomenon which currently 
provides problems to the Standard Model of Galactic CRs. 
The essence of these problems can be easily seen (for a 
detailed discussion see also the papers of other authors 
that cover these topics in the present topical issue). 

The diffusion coefficient and its energy dependence are 
primarily determined by the turbulence spectrum w(k), 
which in the most important cases is given in a power-law 
form w(k) oc k~ m , where k is the wave number. Then one 
has 



D(E) oc E^, fj, = 2-m, 



(8) 



with /i = 1/3, /i = l/2 and \i = for the Kolmogorov, 
Kraichnan and shock wave diffusion, respectively. In the 
Bohm regime with strong turbulence, the diffusion length 
is of the order of the gyroradius and \i = 1. 

It is easy to verify that the space density n p of primary 
protons and nuclei is 



(9) 



It follows from n p (E) ~ Q p (E)T esc (E), where Q P (E) oc 
E~ 1b is the generation (acceleration) rate of primaries, 
i.e. protons and nuclei, and T esc (E) oc D(E)^ 1 is the es- 
cape time from the Galaxy. For non-rclativistic diffusive 
acceleration 7 9 = 2; hence from the index of the ob- 
served spectrum, 7 9 + fi — 2.7, one derives fi = 0.7. At 
very high energy this results in a too large anisotropy, 
5(E) oc D(E) oc and in a too small traversed gram- 
mage, X CT (E) cx l/D(E). 



To solve the problem of small X CI (E) it was sug- 
gested the effect of spallation inside CR sources and re- 
acceleration by second order Fermi mechanism during CR 
propagation [43 1. 

Concerning anisotropy, observations show its constant 
behavior with energy which is in contrast with the simple 
estimate above. This inconsistency is only apparent be- 
cause the result 5(E) oc D(E) refers to an average of the 
anisotropy amplitude computed over many source realiza- 
tions, i.e. over a continuum of sources. In the framework 
of the SM for GCR, the actual value of anisotropy can be 
computed only taking into account the random nature (in 
space and time) of nearby recent SNRs [44j . 

It is easy to estimate the e + /e~ ratio in the framework 
of SM and to compare it with the one measured in the 
PAMELA experiment 0. 

The basic difference with protons and nuclei, considered 
above, is that the lifetime of HE electrons and positrons is 
much shorter than the escape time. It is given by energy 
losses due to synchrotron and inverse Compton radiation, 
T e oc E —1 . For primary HE electrons the rate of production 
(acceleration) is Q e (E) oc E 1g and the space density 



n e (E) ~ Q e (E)r e (E) oc E'^ + ^ 



(10) 



Positrons in the SM GCR arc produced as secondary par- 
ticles in p + p — > 7r + — > e + interactions of CR protons with 
the interstellar gas. Using n p (E) from Eq.©, one finds 
the rate of positron generation as 

Q e +(E) ~ n p (E)a pp n gas c oc E~ { ^+^ . (11) 

The space density of positrons n e + (E) is then obtained as 

n e +(E) ~ Q e +(E)r e oc E~^ + ^ . (12) 
Thus in the framework of the SM of GCR the ratio 

n e +/n e _ oc E^^ (13) 

regularly diminishes with energy increasing. This simple 
estimate is in contrast with the observed increasing of the 
e + /e~ ratio in the PAMELA data [45 1. 



2.3. Standard Model: Problems and Solutions 

The SM of Galactic cosmic rays faces at least two major 
problems. The first one is connected with the diffusion 
mechanism and was already discussed in the previous sec- 
tion l2.2l The second one is connected with the acceleration 
mechanism of CRs and will be discussed below. Here we 
shall address only the contradictions of the SM of GCR 
with observations, demonstrating that in most cases they 
are only apparent. Numerous subdominant phenomena, 
like re-accclcration, which in fact may affect the results in 
a noticeable way, will be omitted. 

As far as diffusion is concerned, the problems arise due 
to the fast increase of the diffusion coefficient with en- 
ergy D(E) oc E^, which results at high energy in a large 
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Figure 1. Fluxes and spectra calculated within the Standard Model |42| for all particles, protons and nuclei. Fluxes are shown as function of 
kinetic energy per nucleon. They are compared with the data of CAPRICE, ATIC-2, JACEE and KASCADE. The position of the knees 
for all nuclei are given by 3Z PeV. The end of Galactic spectrum is given by the iron knee m 80 PeV. At higher energies Galactic 

spectrum becomes very steep. 



anisotropy, 8(E) ex. E 11 , and in a too small traversed gram- 
mage, X cr oc E~^. A possible solution was mentioned in 
section 12.21 

It is interesting that another diffusion problem, the ob- 
servation by PAMELA of the rise with energy in the e + /e~ 
ratio, may be solved within the SM [46| (see, however, the 
criticism of this model in [It}). The basic idea presented 
in [iiH is the production of e + e _ -pairs in pp-collisions in 
situ, i.e. at the shock front. Secondary positrons arc im- 
mediately involved in the process of acceleration, in fact 
more efficiently than primary electrons, because they al- 
ready have a power- law energy spectrum. As a result, the 
ratio e + /e _ rises with energy. Of course, there were pro- 
posed other mechanisms of positron production, e.g. from 
pulsars 0, [49[, from Dark Matter annihilation etc. (see 
[501 ] for review and references), but it is important that 
the positron excess can be successfully explained in the 
framework of the SM GCR. 

The difference in spectrum slopes for protons and He- 
lium was considered for a long time as a serious prob- 
lem for shock acceleration, which is the most important 
component of the SM. This difference is presented in the 
most precise way by the PAMELA experiment 51] ; the 
spectrum of Helium is harder. The slope difference is 



small, but statistically it is reliably provided, 7 P — 7h c = 
0.1017 ±0.0014. 

Theoretically, the shock-wave acceleration does not dis- 
tinguish among nuclei with equal rigidities. The predicted 
spectrum is a decreasing power-law with an exponent de- 
termined by the shock properties, i.e. the Mach number. 
However, according to PAMELA, spectra are not exactly 
power-law and the effective exponents are also different: 
2.8 and 2.7 for protons and Helium, respectively The un- 
derstanding of these phenomena needs the inclusion of par- 
ticle injection in the acceleration regime and the account 
of the escaping particles from the accelerator. We shall 
follow here the paper by Drury [52| with a clear physical 
analysis and the works 5J3, l54| for calculations. 

In the standard approach the spectrum is calculated for 
particles accompanying the shock front, i.e. located within 
the acceleration length upstream and downstream. This 
spectrum differs from that of the escaping and captured 
particles. The latter drift downstream and are trapped 
there until the shock dies. They suffer adiabatic energy 
losses and do not contribute to the HE part of the ob- 
served spectrum. For escaping particles the exit model 
must be specified. In section 12.11 we considered the exit 
mechanism from (37| : at each time t particles accelerated 
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up to maximum energy £ , max (i) freely escape from the 
shell. The total spectrum of escaping particles may be 
calculated by an integration over time. The highest en- 
ergy particles escape at the beginning of the Sedov phase; 
the lower energy ones exit progressively later. The chem- 
ical composition changes with time because swept-up gas 
changes its chemical composition and ionization state. It 
means that the energy spectrum of nuclei with fixed (Z, 
A) is characterized by its own j g . 

According to section 12.11 injection plays a similar role. 
An account for propagation through interstellar medium 
can also result in a difference in the spectral indices of 
protons and Helium. 
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2. 4- Beyond the SM: New component of the Galactic Cos- 
mic Rays? 

The observed knees are interpreted as the end of the 
Galactic CRs. Positions of knees define maximum acceler- 
ation energies, different for protons and nuclei, see Eq. (gj). 
Though this picture looks now very natural, the way to its 
adoption was rather dramatic. Soon after the discovery of 
DSA :34| it was understood (HBJ that the maximum acceler- 
ation energy in this approach does not exceed Z x 10 TeV, 
much below the knee. Only later after the works by Bell 
341 ] . and Bell and Lucek |35j, it was realized that an in- 
creasing of the magnetic field due to the Weibel instability 
may enhance the maximum energy up to the needed level 
£™<™ w o.l EeV given by Eq. (gj). This story teaches us 
to be careful with assumptions about higher £ max . Nev- 
ertheless, reasonable mechanisms to enhance E max have 
been recently found. 

The motivation for an additional component of the 
GCRs with -Emax up to 10 EeV is given by the interpre- 
tation of the ankle observed according to the latest data 
at (3 - 4) EeV, as the transition from Galactic to extra- 
galactic CRs at E tI ~ (3 — 10) EeV. It was proposed in 
IH, Hi], El] that such an additional component may appear 
due to some unspecified sources different from SNR. 

A plausible idea for increasing E m!iX is given by the ex- 
plosion of SN in the stellar wind left behind by presuper- 
nova [Hil, [6(| ■ The large magnetic field in the wind and 
the presence of heavy nuclei result in a larger £ max . In 
particular in [59| this effect is maximized by considering 
the Wolf-Rayet stars as presupcrnovae. 



In a more general way this idea was considered in [61 1. 
In this work an approximate relation between i? m ax and 
a SN total energy output £ is obtained as logEmax ~ 
0.52 log £51 + 2.1 TeV for the protons, where £51 is the 
total energy output in units 10 51 erg. The distribution of 
SNc over £ max is obtained using the distribution of SNe 
over £51. The high energy tail of the distribution gives 
sources with enhanced E max , and thus the diffuse spec- 
trum at large E extends up to E ~ 1 EeV for Iron nuclei. 

In [62] ] a similar approach was used for the calculation of 
GCR diffuse flux from SNRs with large E maK . SNRs were 



Figure 2. Diffuse spectra of Galactic CRs produced by four groups 
of SNRs, labeled as SN la, SN IIP, SN Ib/c and SN lib (see [gj). 
Along the ordinate axis the plotted quantity is p 2 n(p)/ E$n, where p 
is the particle momentum, n(p) is proportional to the total number of 
accelerated particles with momentum p produced by SN and -Bsn is 
the total energy of SN explosion. Note that we generalized the orig- 
inal plot from 1621 introducing dimcnsionlcss rigidity R = p/(Zm^) 
instead of dimcnsionlcss proton momentum p/mjy, as the authors 
recommended to do for nuclei. As seen from the plot, in the maxi- 
mum of the curve the rigidity reaches i? m ax — 3 x 10 7 , i.e. 0.78 EeV 
for Iron (see text). 



subdivided into four groups according to the type of ex- 
plosion: SNR la, SNR Ilb/c, SNR IIP and SNR lib. Each 
group of SNRs was characterized by fixed parameters of 
SN explosion and presupernova stellar wind, in which the 
shock propagates, accelerating particles. In Fig. [2] spec- 
tra of accelerated particles and i? max arc shown for each 
calculated group. The group SN lib represents a new HE 
component with a higher £ max . This energy corresponds 
to the shock propagation through a presupernova stellar 
wind. For the evolution of SN lib we use calculations pub- 
lished in [gl] . The presupernova of SN lib is a red super- 
giant with mass 18Mq on the Main Sequence. The stellar 
wind is presented by Helium [Z = 2) and at a smaller 
extent by Oxygen (Z — 8), since metals have very small 
abundances in presupernova stellar wind. The metals pro- 
duced by SN explosion work as a piston, and the shock 
propagates ahead of this layer. Using Oxygen we obtain at 
the maximum of the curve -Emax = ZR max w 0.2 EeV, and 
for Iron which has very low abundance i? max « 0.8 EeV, 
still below the ankle. 

The next logical step in the discussion of SN explosions 
as sources of a new HE Galactic component is given by 
the hypernova explosion, a very rare (one per 10 5 years 
in a galaxy) but very powerful event which can alone fill 
the Milky Way with cosmic rays of the observed density. 
Hypernovae are assumed to be the sources of GRBs. Such 
a model for GCR was considered in [@4j|. The explosion 
occurred 2 x 10 5 yr ago at distance about 500 pc from the 
Sun. At present the highest energy particles passed the 
Sun, while the lower energy ones are still confined there. 
In this model the transition to extragalactic CRs occurs 
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at E ~ 0.2 EcV, too low for the "ankle transition". 

An extreme case is considered in [65[. A sample of 1000 
hypernovae are exploding randomly in the Milky Way with 
time interval 10 5 years. They produce cosmic rays with a 
mass composition of 10% protons and 90% Iron nuclei with 
maximum energy 100 EeV. Protons, propagating rectilin- 
early, leave the Galaxy starting from E ~ 1 EeV, while 
Iron nuclei diffuse in the Galactic magnetic holds and well 
explain the energy spectrum of the PAO starting from 
10 EcV. At lower energy the spectrum and mass composi- 
tion are explained as a mixture of protons and Iron nuclei. 
This is the case when hypernovae can explain the whole 
spectrum of UHECR at E > 2 EcV by the Galactic CR 
component. 

To conclude, the ankle model needs a new high energy 
Galactic component at 0.1 EeV < E < 10 EeV with a 
heavy mass composition. Though in principle it looks pos- 
sible, for SNe such a case was not found. Probably other 
Galactic sources, such as pulsars [3l|, could fit this hope 
better. 

3. Extragalactic UHECR 

In this section the study of the transition will be pro- 
ceeded moving from the end of the UHECR spectrum 
(from GZK cutoff in the case of protons) towards lower 
energies. 

In the case of protons there are two spectral signa- 
tures of their propagation through CMB: pair-production 
dip, which is a rather faint feature caused by the pair- 
production process given by Eq. (fTJ), and a sharp steep- 
ening of the spectrum caused by pion photo-production 
process 

P + 7cmb -> 7T + X, (14) 

called GZK cutoff [J Q. The GZK cutoff position is 
roughly defined by the energy where the pair-production 
energy loss in Eq. ((TJ) and the pion production one of 
become equal, namely at Eqzk 



12 



Eq. ([14]) become equal, namely at E'gzk — 50 EeV 

In the case of UHE nuclei the situation is different. For 
heavy enough nuclei A > 10, taking into account only the 
CMB background, the steepening of the spectrum starts 
at the energy where adiabatic energy losses (due to the 
expansion of the universe) and losses due to pair produc- 
tion become equal (see energy losses in 0), namely at 
E r-j 20 — 50 EcV. The pair production process experi- 
enced by nuclei is affected only by the CMB radiation, 
while the photo-disintegration process takes place also on 
the Extragalactic Background Light (EBL) [67[. The lat- 
ter is the process in which a nucleus with atomic mass 
number A loses one or more nucleons interacting with the 
CMB or EBL: 



A + jcmb.ebl ->(A- nN) + nN 



(15) 
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Figure 3. Modification factors for different generation indices f g . 

gether with the pair production process, produces a steep- 
ening in the observed spectrum. The exact position of 
the flux suppression and its shape depend on the nuclei 
species as well as on the details of the cosmological evolu- 
tion of the EBL which, in contrast to CMB, is not known 
analytically, being model dependent (67}. In general the 
steepening of the nuclei spectrum is not as sharp as the 
GZK and it occurs at lower energies (66|. 

Thus, the GZK cutoff is the most spectacular feature 
in the UHECR spectrum. However, it is valid only for 
UHE protons and the predicted shape of the GZK fea- 
ture is model-dependent: it depends on a possible lo- 
cal over-density or deficit of sources and can be mim- 
icked by low values of the maximum energy E^ x , that 
sources can provide, In contrast, the pair-production dip 
is practically model-independent mainly because protons 
contributing to this energy range arrive from distant cos- 
mological sources. 

3.1. Pair-production dip 

Being a quite faint feature, the pair-production dip is 
not seen well in the spectrum presented as log J(E) versus 
log_E. It is more pronounced when analyzed in terms of 
the modification factor 



V (E) = J p (E)/J»™(E), 



(16) 



The dominant process is one nucleon emission (n = 1) as 
discussed in [661 ] - The photo-disintegration of nuclei, to- 



which is defined as the ratio of the spectrum J p (E), cal- 
culated with all energy losses taken into account, and the 
unmodified spectrum J p nm (E), where only adiabatic en- 
ergy losses are included. 

Theoretical modification factors calculated for different 
source generation indices 7 S are presented in Fig. [3] If one 
includes in the calculation of J P {E) only adiabatic energy 
losses, then, according to its definition, Tj(E) = 1 (dash- 
dot line). When the e + e~-production energy loss is ad- 
ditionally included, one obtains the rj(E) spectrum shown 
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Figure 4. The predicted pair-production dip in comparison with Akeno-AGASA, HiRes, Yakutsk and Telescope Array (MD and SD) data 
[17|I . All these experiments confirm the dip behavior with good accuracy, including also the data of Fly 's Eye 11711 (not presented here) . 



in Fig. by the curves labeled "e + e~". With the pion 
photo-production process being also included, the GZK 
feature (curves "total") appears. The observable part of 
the dip extends from the beginning of the GZK cutoff at 
E w 40 EcV down to E « 1 EcV, where 77 w 1. It has 
two fattcnings: one at energy E^ ~ 10 EeV and the other 
at Ef, ~ 1 EcV. The former automatically produces the 
ankle (see Fig. [4]) and the latter provides an intersection 
of the flat extragalactic spectrum at E < 1 EcV with the 
steeper Galactic one. Since many non-essential factors in 
the numerator and denominator of Eq. (|16p compensate 
or cancel each other, the dip in terms of the modification 
factor is a less model dependent physical quantity than the 
spectrum. In fact the dip is determined mostly by the in- 
teraction of protons with the CMB photons and it depends 
mainly on the CMB spectrum and the differential cross- 
section of e + e~ pair-production. In particular it depends 
weakly on the spectral index of the generation spectrum. 
In Fig. [3] curves are plotted for 2.1 < j g < 3.0 with a step 
A7 ff = 0.1, and uncertainties are seen as thickness of the 
curves. 

Modification factors in Fig. [3] are presented for the case 
without cosmological evolution of the sources, which is 
usually described in the injection spectrum by a factor 
(1 + z) m up to z max . The inclusion of evolution may no- 
ticeably change the modification factor, but in fact it al- 



lows to improve the agreement of the dip with data due to 
the additional free parameters m and z max (see Fig. 14 of 
Ref. [H). 

Thus, a remarkable property of the dip in terms of mod- 
ification factor is its universality. The dimensionless quan- 
tity rj(E) remains the same with various physical phenom- 
ena being included in calculations 15|, [l6( : discreteness in 



the source distribution (distance between sources may vary 
from 1 Mpc to 60 Mpc), different modes of propagation 
(from rectilinear to diffusive), local overdensity or deficit of 
sources, large-scale inhomogeneities in the sources distri- 
bution, some regimes of cosmological evolution of sources 
(most notably those observed for AGN) and fluctuations 
in the interactions. The phenomenon which modifies the 
dip significantly is the possible presence of more than 15% 
of nuclei in the primary radiation. Therefore, the shape 
of the proton dip in terms of modification factor is deter- 
mined mostly by the interaction with the CMB. 

Above the theoretical modification factor was discussed. 
The observed modification factor, according to definition, 
is given by the ratio of the observed flux J bs(E) and un- 
modified spectrum J un m(E) cx E 1-79 , defined up to nor- 
malization as: 

Vohs<x J ohs {E)/E~^. (17) 
Here 7 g is the exponent of the generation spectrum 
Qgcn(Eg) cx Eg 7s in terms of initial proton energies E g . 
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Figure 5. Mass composition and GZK cutoff as measured by the HiRes detector. In two upper panels (X max ) (left) and RMS (right) are 
presented as function of the energy. Both agree with a pure proton composition, shown by curves labeled 'proton'. The left- lower panel 
shows differential energy spectrum in terms of the modification factor. One can see a good agreement with the predicted shape of the GZK 
cutoff. The right-lower panel shows the quantity E\/2 in the integral spectrum. This energy, a characteristic of the GZK cutoff, is found as 
^1/2 = 10 19 ' 73±0 ' 07 eV in good agreement with theoretical prediction E 1 / 2 = 10 19 ' 72 eV (see the text). 



Fig.rj]shows that both the pair production dip and the be- 
ginning of the GZK cutoff up to 80 EeV are well confirmed 
by experimental data of Akeno-AGASA, HiRes, Yakutsk 
and TA [l7|. 

The comparison of the theoretical dip with observa- 
tional data includes only two free parameters: exponent 
of the power-law generation spectrum j g (the best fit cor- 
responds to 7 g = 2.6 — 2.7) and the normalization con- 
stant to fit the e + e~-production dip to the measured flux. 
The number of energy bins in the different experiments is 
20 - 22. The fit is characterized by x 2 /d.o.f. = 1.0-1.2 
for AGASA, HiRes and Yakutsk data. For this fit we used 
the modification factor without cosmological evolution of 
sources. As was explained above, using two additional evo- 
lutionary parameters, m and z ma x, one can further improve 
the agreement. 

In Fig. [5]one can see that at E < 0.6 EeV the exper- 
imental modification factor, as measured by Akeno and 
HiRes, exceeds the theoretical modification factor. Since 
by definition the modification factor must be less than one, 
this excess signals the appearance of a new component of 
cosmic rays at E < E tr ~ 0.6 EeV, which can be noth- 
ing else but the Galactic cosmic rays. This interpretation 
is confirmed by the transition of a heavy component to 
protons in the upper-left panel, that with good accuracy 
occurred at the same energy. Thus, according to HiRes 



data the transition from extragalactic to Galactic cosmic 
rays, occurs at energy E ir ~ 0.6 EeV and is accomplished 
at E ~ Eb rj 1 EeV (see left-upper panel in Fig. [5] as 
example). 

3.2. GZK cutoff. 

The two largest Extensive Air Shower (EAS) detectors, 
HiRes and Pierre Auger Observatory 0, Q , have observed 
a sharp steepening in the UHECR spectrum at E > (30 — 
50) EeV. Both collaborations claimed that the observed 
steepening is consistent with the GZK cutoff. But as a 
matter of fact, there is a dramatic conflict between these 
two results, which still leaves the problem open. 

In this subsection we analyze the data of HiRes which 
provide a strong evidence in favor of the GZK cutoff. 
These data are supported also by the TA observations 0| ■ 
The data of PAO will be considered in the next subsection. 

To interpret convincingly the spectrum steepening as 
the GZK cutoff one must prove that: 

1. energy scale of the cutoff and its shape correspond to 
theoretical predictions, 

2. the measured chemical composition is strongly domi- 
nated by protons. 

Mass composition of the CR spectrum in turn may be 
measured with the help of 
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Figure 6. PAO spectrum I68t l69l| (filled squares) compared with HiRes/TA spectra and with calculated spectrum (full line) with parameters 
indicated in the figures. In the left panel the spectra are shown with energy scales as determined in the experiments. In the right panel the 
data of all three experiments are recalibrated using the pair-production dip as a standard candle. For HiRes and TA the recalibration factor 
is 1.0, since these detectors describe well the dip. The PAO recalibration factor should be 1.22 to fit the dip shape. It is not trivial that after 
this recalibration the PAO flux coincides well with HiRcs and TA fluxes. However, the contradiction of three PAO bins at 40 — 50 EeV with 
the predicted GZK shape remains. 




E[eV] 





QGSJET01 


■ Auger 09 




QGSJETII 


proton 




< tft m; ( 








'•M : 












iron 



10* 



10" 



E[eV] 



Figure 7. Auger data [70] for (X max) as function of the energy (left panel) and for RMS(-Xmax), the width o f th e distribution over X max , 
(right panel). The calculated values for protons and Iron are given according to QGSJET01 I7lll and QGSJETII [72f] models. One can see from 
the right panel that RMS distribution becomes more narrow with increasing energy which implies a progressively heavier mass composition. 



(I max ) (E), the average depth of atmosphere in 
g/cm 2 , where a shower with energy E reaches its max- 
imum, and 



• RMS(V max ), 
over X max . 



which is the width of the distribution 



These values measured by HiRes are displayed in Fig. [5] 
From the left-upper panel of Fig. [5] one can see that the 
chemical composition changes from very heavy elements, 
probably Iron, at E ~ 0.1 EeV, (data of HiRes-MIA Jz3] ) 
to protons at E ~ 1 EeV (data of HiRes [I?} ) ■ Such evo- 
lution fully corresponds to the SM for GCR (see Fig. [1]). 
RMS(V max ), a very sensitive tool for mass composition, 
also provides evidence for a proton-dominated composi- 
tion at E > 1 EeV and up to the highest energies (see 
upper-right panel of Fig. [5]). Differential energy spectrum 
of the GZK feature in the modification factor form (left- 
lower panel) is in a reasonably good agreement with the 
theoretical prediction, though better statistics at higher 
energies is still needed for a final conclusion. 

The integral energy spectrum of UHE protons, J p (> E), 
her specific i 

h/2 m it 

calculated integral spectrum below 50 EeV is well approx- 



has another specific characteristic of the GZK cutoff, the 
energy E 1 / 2 UM- It is based on the observation that the 



imatcd by a power-law function: J p (> E) oc E~^. At high 
energy this spectrum steepens due to the GZK effect. The 
energy where this steep part of the spectrum equals to half 
of its power-law extrapolation, J p (> E) = KE^^ , defines 
the value of Ei/ 2 . This quantity is found to be practically 
model-independent; it equals to Ei/ 2 = 10 19 72 eV ~ 52.5 
EeV [12j. Fig. [5] demonstrates how the HiRes collabo- 
ration found Ex/ 2 from observational data [74j]. The ra- 
tio of the measured integral spectrum J(> E) and the 
low-energy power-law approximation KE~ 7 was plotted 
as a function of energy. This ratio is practically constant 
in the energy interval 0.3 — 40 EeV, indicating that the 
power-law approximation is a good fit, indeed. At higher 
energy the ratio falls down and intersects the horizon- 
tal line 0.5 at the energy defined as Ei/ 2 . It results in 



E 



1/2 



10 



19.73±0.07 



eV, in excellent agreement with the 



predicted value. 



Thus, one may conclude that the HiRcs data presented 
in Fig. [5] indicate the proton-dominated chemical composi- 
tion and the presence of the GZK cutoff in both differential 
and integral spectra. The conclusion about proton com- 
position is further supported by the recent TA data 17 [. 
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Figure 8. Left panel: Energy spectrum in the two-component model with protons and Iron nuclei for homogeneous distribution of sources and 
with 7 9 = 2.0 and -B ma x = 4 EeV. Right panel: As in left panel but with a diffusion cutoff. The gap is expected to be filled by intermediate 
mass nuclei. 



3.3. PAO data: energy spectrum and mass composition 

The spectrum of the Pierre Auger Observatory, shown 
by filled squares in the left panel of Fig. |6l differs from the 
HiRcs and Telescope Array spectra both in absolute nor- 
malization and in shape. However, one is allowed to shift 
the energy scale within systematic energy errors. Since the 
energy position of the pair-production dip is ri gidl y fixed 
by the interaction with CMB, in Ref. [Iff and [16J it was 
proposed to use the dip as an energy calibrator ('standard 
candle'). The data of HiRes and TA have recalibration 
factors fc ca i = 1.0, since these data agree well with the 
energy position and shape of the pair-production dip (see 
left panel of Fig. [5]). The Auger data (left panel of Fig. \6§ 
disagree with the dip, but with the recalibration factor 
fc ca i ~ 1.22 a reasonable agreement can be reached. It is 
far non-trivial that simultaneously the Auger flux agrees 
well with the fluxes of HiRes and TA. This procedure was 
applied in [IB], [l6| imposing the condition of minimizing 
X 2 after recalibration of the data of AGASA, HiRes and 
Yakutsk detectors. A similar procedure based on the dip 
energy calibration was also used by Blumer et al. [IlJ for all 
experiments with recalibration factors 1.0, 1.2, 0.75, and 
0.625 for the HiRes, PAO, AGASA and Yakutsk detectors, 
respectively; a good agreement for all fluxes was demon- 
strated. This procedure was also successfully applied in 
75]. 

The apparent coincidence of the PAO and HiRcs/TA 
spectra is related to the low energy part of the energy 
spectrum in Fig. [51 At higher energies statistical uncer- 
tainties are too large to distinguish between the spectra. 

While the HiRes and TA spectra are compatible with the 
GZK cutoff, the Auger spectrum is not. The steepening 
in the upscaled PAO spectrum starts at energy E < 40 
EeV, lower than Eqzk — 50 EeV, and in three successive 
energy bins the PAO flux is significantly lower than one 
predicted for the GZK shape as shown in the right panel of 
Fig. |51 We could not reconcile the PAO cutoff shape with 
the GZK behavior by including in calculations different 
generation indices j g , evolution regimes, low acceleration 
maximum energy -E max , local over-density of sources etc. 

This disagreement is quite natural taking into account 
the mass composition measured by PAO (see Fig. [7]), 
which dramatically differs from that of HiRes shown in 



Fig. [5] The dependence of the RMS(V max ) on energy, 
which steadily decreases and approaches the Iron value at 
E ~ 35 EeV, is a strong example of this disagreement. 
Low RMS, i.e. small fluctuations, is a typical and reliable 
feature of a heavy nuclei composition. These data are fur- 
ther strengthened by other PAO measurements provided 
by surface detectors. They allow to extract two mass- 
composition dependent quantities: the atmospheric depth 
(-X"max)) where muon-production rate reaches maximum, 
and maximum zenith angle # max determined by the sig- 
nal rise-time in surface Cerenkov detectors. Measurements 
of both quantities confirm the heavy mass composition 
and its dependence on energy obtained with the help of 
(A max ) (E) and RMS(A max ) and (z^. The soon ex- 
pected data on muon flux from the Auger Muon and Infill 
Ground Array (AMIGA) [zl| will further clarify the mass 
composition. 

An early steepening of the spectrum observed by the 
PAO can be explained as a property of nuclei propagation. 
However, it is a problem to explain simultaneously both 
spectrum and mass composition of PAO. 

A wide class of mixed composition models pjj fails to 
fit the PAO mass composition at the highest energies. 
Most of them assume a source mass composition close to 
Galactic CRs, i.e. enhanced by protons. The observed 
mass composition in these models becomes heavier only 
at E > 50 EeV, when protons disappear due to the GZK 
cutoff [ii]. In the paper {HI, by Allard ct al. (2008), an 



exceptional case of primaries being pure Iron nuclei was 
considered, though authors refer to this possibility as to 
an unnatural one. This model (see also [2j fig. 4 and [80j |) 
explains the observed Auger energy spectrum. 

Mass composition becoming progressively heavier with 
increasing energy appears in the rigidity-dependent accel- 
eration, when the maximum attainable energy for a nu- 
cleus with charge number Zi is £™ ax = Zi_E™ ax and the 
contribution of nuclei with charge Z' < Z to this energy 
disappears, while heavier nuclei with larger Z exist. In 



811 was built. 



this way the 'disappointing model' 

The disappointing model is based on the following ob- 
servations and assumptions: 

• According to the HiRes (Fig. [5} and PAO (Fig. [7J 
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data, the observed primaries at energy (1 — 3) EeV 
are predominantly protons. 

• It is assumed that they have extragalactic origin. 

• The particles at higher energies are extragalactic nu- 
clei with charge number Z steadily increasing with 
energy. 

• Acceleration of CRs in sources is rigidity dependent, 



E? 



= Zi x £™ ax . 



It was demonstrated in [81j that to avoid a proton dom- 
inance at the highest energies, one must assume that the 
maximum energy of the accelerated protons is limited, 
£™ ax < (4 — 10) EeV. This conclusion is valid for a large 
range of generation indices, 7 9 ~ 2.0 — 2.8, and for a 
wide range of cosmological evolution parameters. The cal- 
culated proton and nuclei energy spectra for j g = 2.0, 
£ max = Z x 4 EeV and without cosmological evolution are 
shown in Fig. [5] In the left panel the two component model 
(protons and Iron) is presented. In the right panel of Fig. [5] 
intermediate primary nuclei are included in the framework 
of the diffusive propagation through intergalactic magnetic 
fields (see [8l| for details). The rigidity dependent accel- 
eration, £'™ ax — Zi x E"™ ax , provides a steadily increasing 
charge number Z as energy increases. 

The disappointing model as it is presented here and in 
8l| is not complete. The energy spectra must be cal- 
culated in diffusive approximation for different distances 
between sources. Energy dependence of (V max ) (E) and 
RMS(V max ) must be consistent with calculated energy 
spectrum. 

This model was called 'disappointing ' because of lack of 
many signatures predicted in proton-dominated models, 
such as cosmogenic neutrino production and correlation of 
CR arrival directions with distant sources. 

Recently a similar model [82[, based on the rigidity- 
dependent acceleration, was proposed for making the 
chemical composition heavier with energy increasing. This 
model uses the reasonable assumption [8^, [l6| that the 
space density n s of UHECR sources decreases with in- 
creasing of the maximum acceleration energy. The idea 
may be explained in the following way. 

Consider accelerators which produce protons (p) and nu- 
clei (i) with maximum energy £™ ax and E™ &x = Zi x 
£™ ax , respectively. The space density of sources is as- 
sumed to depend on E max , 



n s (E p ) CX (E p ) 



max ^ 
P 



-0 



Generation rate of protons with energy E per unit space 
volume is proportional to 



Q p (£)oc [ dE p ™*(E™*)- fj oc K 
Je 



(18) 



because sources with i?" lax < E do not contribute. For 
nuclei i 



Qi(E) oc / dE^m 

J E/Zi 



(19) 



(sources with £™ ax < E/Z, 
contribute). Then from Eqs. 



Zrmax 

and iff 



i(E)/Q p (E) oc Z- 



(9—1 



< E, do not 



(20) 



Eq. (|20|) shows that the production of nuclei is enhanced 
by a factor Z®~ 1 ; it can affect the mass composition cal- 
culations in the mixed composition models. 

4. Three transition models 

In this section we discuss three models of transition from 
Galactic to extragalactic CRs: ankle, dip and mixed com- 
position models. One feature is common for all of them: 
the transition is described as an intersection of a steep 
Galactic spectrum with a more flat extragalactic one. The 
agreement with the SM for GCR (see subsection l2.1[) is one 
more criterion which these models have to respect. Ac- 
cording to the Standard Model, the benchmark of the end 
of GCR is given by the Iron knee at energy Ep e rj 80 PeV 
(see Fig. [TJ and at E > Ep e it has an exponential cutoff. 

Motivated by the interpretation of the ankle as the tran- 
sition to extragalactic CR at E^ ~ (3 — 10) EeV, one has 
to assume [Hfl [H3, 58 1 an additional component of GCR 
accelerated up to energies much beyond the Iron knee (see 
subsection 12.41) . 

Observational data which have the power either to con- 
firm or to reject each transition model include energy 
spectrum, elongation curve (V max ) (E), RMS(V max ) and 
anisotropy. Below these models will be discussed in the 
historical succession of appearance: ankle, dip and mixed 
composition. 

4--1- Ankle model 

This is the traditional model based on the interpreta- 
tion of the ankle as the spectrum feature where transi- 
tion occurs (see e.g. [22 



23j). In fact, this is a very 



natural model since transition occurs because the extra- 
galactic component is very hard. This component is as- 
sumed to have a pure proton composition with a flat gen- 
eration spectrum Q ox t r . P oc E~ 2 valid for non-relativistic 
shock acceleration. Energy losses modify the spectrum in- 
significantly at E < 40 EeV. The beginning of the ankle 
at E° hs ~ (3 — 4) EeV corresponds to the energy where 
fluxes of Galactic and extragalactic CRs get equal. Thus, 
the Galactic CRs should be presented by an additional 
component accelerated up to energy at least by factor 
30 — 40 times higher than the maximum energy in the 
Standard Model. In the majority of the ankle models, e.g. 
56, 57, 8f| 23|, the large fraction of the observed cosmic 
rays has a Galactic origin at E > 10 EeV. To facilitate 
the acceleration problem one should assume a heavy-nuclei 
composition of the new component. 

The transition at the ankle is illustrated by the right 
panel of Fig. [9l The curve "extr.p" presents the calcu- 
lated extragalactic flux of protons and the dash-dot line 
gives the Galactic CR spectrum. The latter is obtained 
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Figure 9. Transition in the dip (left panel) and ankle (right panel) models. In both cases solid lines give the calculated extragalactic proton 
spectrum and dashed lines describe the Galactic iron spectrum. Etr is the energy of intersection of the Galactic and extragalactic spectra 
and Ep e gives the position of the iron knee. E^ rs 1 EeV in the left panel is the energy where the transition from Galactic to extragalactic 
CRs is completed. 



by subtracting the extragalactic component from the to- 
tal observed flux following the procedure first suggested in 
[8?] ]. The observed dip in the spectrum may be explained 
by the Hill-Schramm's mechanism (lcj |. 

Another problem of the ankle model is the contradic- 
tion with the measured average depth of EAS maximum, 
(X max ) (E), in the energy range (1 — 5) EeV. While all 
data, including both HiRcs and PAO, show proton or light 
nuclei composition here, the ankle model needs a heavy 
Galactic component, predicting too small (A" max ) (E) in 
contradiction with observations (see the right panel of 
Fig. [10] and right panel of Fig. 4 in HH). This contra- 
diction is found also in 
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The ankle at an energy higher than 3 EeV contradicts 
also the anisotropy calculated in [89| . 



4-2. Dip model 

The dip model is based on the assumption that UHECR 
at E > 1 EeV arc mostly extragalactic protons. This 
assumption is confirmed by the observation of the pair- 
production dip in the energy range (1 — 40) EeV and the 
beginning of the GZK cutoff in the energy range (40 — 
100) EeV (see Fig. [5]). Both features are signatures of 
a proton dominated spectrum. As was discussed above 
the shape of the dip allows an admixture of light nuclei, 
though not more than 15%. The transition from Galactic 



to extragalactic CRs occurs as the intersection at E tl - ~ 

0. 5 EeV of the steep Galactic component (dashed line in 
the left panel of Fig. ^ with the flat extragalactic proton 
component shown by 'extr.p' curve (this curve appears to 
be falling down since the spectrum is multiplied by E 25 ). 
The flatness of the extragalactic spectrum is provided by 
the distribution of sources over the maximum acceleration 
energy n s (E niax ) oc E~@ x (see below) or, in the case of 
diffusive propagation, by the 'magnetic flattening' li[ 9(| 
9l| . The transition is completed at energy Ef, « 1 EeV, 

1. e. it occurs at a visible feature in the CR spectrum known 
as the 'second knee', at energy between E tI and Ef, in the 
left panel of Fig. M 



The basic features of the dip model wee as follows [14, 



• The primary flux is strongly proton-dominated. 

• Sources are e.g. AGN 1921 with a neutron mechanism 
for particle escape 93j, l94| which provides a pure pro- 
ton generation spectrum. Another case of enhance- 
ment in the proton fraction at relativistic shock ac- 
celeration is given in (l6j . 

• To reproduce the observed shape of the dip, indices 
for the generation rate per unit co-moving volume 
Q(E) oc E~ lg in models without evolution are to be 
j g = 2.6 — 2.7. Such steep spectra can be obtained for 
the case of usual injection q(E) oc (E~ 2 — E~ 2 - 3 ) from 
a single source assuming sources distribution over 
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Figure 10. The calculated elongation curves (X max ) (E) for the dip model (two left panels) and ankle model (right panels). The calculated 
curves (X max ) (E) are shown by the thick solid lines for QGSJET01 [fjl model of interaction, by the thick dashed lines for QGSJET-II [72I] , 
and by dotted lines for SIBYLL [gj. The data points are measurements of HiRes-Mia (filled triangles), HiRes and PAO (both filled boxes). 
The PAO data in upper panels with systematic errors, shown by the thin curves, are taken from |85| . The lines 'p' and 'Fe' present the 
elongation curves for proton and Iron which are used for calculations of the model elongation curves in each panel. For PAO data they are 
theoretical curves, for HiRes they include cuts, detector's properties etc. and are taken from curves 'iron' and 'proton' in the upper-left panel 
of Fig. [5] As main result of these plots one may notice the great discrepancy of the ankle model with the data. 



maximum acceleration energies n s (E max ) oc E m g K 



83j, |16[, where n s is the space density of sources. 



• The generation index j g — 2.6 — 2.7 is the main fac- 
tor which provides the agreement of proton spectrum 
with the pair-production dip and the GZK cutoff (see 
subsection 13. ip . The cosmological evolution of the 
sources can be easily included, affecting mostly the 
low energy part of the spectrum. Inclusion of addi- 
tional parameters allows to improve the fit. In partic- 
ular, in [l5| the spectrum was calculated with an ac- 
count for cosmological evolution of AGN as it follows 
from X-ray observations. The calculated spectrum 
shown in Fig. 14 of Ref. [l5| is in excellent agreement 
with data. 

• The mode of propagation (from rectilinear to diffu- 
sive) is important in two energy ranges: at low ener- 
gies E < 1 EeV and at E > 60 EcV. In both regimes 
the distance between sources becomes important in 



the diffusive regime. In the low-energy range it makes 
the spectrum at E < E tr more flat, "magnetic flat- 
tening" shown by curve 'extr.p' in Fig. [5J And at 
E > 60 EeV the GZK shape becomes more steep in 
the case of large distance between sources (see [95| - 



96|). 



The confirmation of the dip model follows from: (i) the 
agreement of the dip energy spectrum shape (see Fig. U 
and Fig. with observations, (ii) the agreement of fluxes 
of all experiments after the dip-based energ y re calibration 
(see Fig.Uand Refs. pjl, Bliimer et al. from [13, [z^), (m) 
the agreement of (X max ) (E) with the bulk of observational 
data in the left panel of Fig. [TTJ1 excluding only the highest 
energy PAO experimental points. 

A characteristic feature of the dip model is the sharp 
transition from Galactic Iron to extragalactic protons. It 
results in a steep increase of the (A max ) (E) with energy 
at E below 1 EeV (see left panels of Fig. [TO]). This feature 
differentiates the dip model from the ankle model, where 
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the increase of (X max ) (E) is less steep (see the right panel 
of Fig. [TU]) . The difference is due to less hard proton spec- 
trum at E < 1 EeV in the ankle model, see Fig. |H1 

The key observation to accept or reject the dip model is 
the chemical composition of UHECR at E > 1 EeV. In the 
case of a substantial admixture of nuclei in the spectrum 
(> 20%) the dip model should be rejected. 

At present the dip model is confirmed by the data of 
HiRes by strong proton-dominance at E > 1 EeV, and is in 
contradiction with the Auger measurements of (A max ) (E) 
and RMS at E > 4 EeV. 

4-3. Mixed composition model 

The main concept of the mixed composition model (see 
Allard et al. [28| , 97 1 , [98| ) is based on the argument that 



any acceleration mechanism operating in gas involves dif- 
ferent nuclei in acceleration process and thus the primary 
flux must have a mixed composition. 

The basic features of the mixed composition model are 
as follows [II [§3, [13: 



• It is assumed that the Galactic component of CRs 
above 0.1 EeV has a pure iron composition while 
the extragalactic CRs have a mixed-nuclei composi- 
tion. Transition can occur either at ankle or below it, 
down to the second knee. Its position is parameter- 
determined. 

• In its basic variant [28| . the source composition of 
extragalactic CRs is assumed to be almost the same 
as for Galactic CRs, with protons and Helium being 
the dominant components. However, a pure species 
composition is also possible in the framework of this 
model, in particular pure Iron, for understanding the 
physical effects. 

• The source energy spectra are taken as power-law 
Q(E) oc E~ la with a generation index 7 9 = 2.1 — 2.3. 
The maximum acceleration energy is assumed to be 
rigidity dependent with a maximum energy for nu- 
clei i being E™- ax = ZiE™ ax , where Zi is the nucleus 
atomic number. The maximum energy of protons is 
typically taken E 7 p nax = 300 EeV. 

• The cosmological evolution of sources is described by 
the factor k cv = (1 + z) m for a wide range of regimes: 
m = in the absence of evolution, m = 3 like for star 
formation evolution and m = 4 for a strong cosmolog- 
ical evolution. These regimes are typically considered 
only for small redshifts z < 1 or 1.3 and for larger z 
the frozen evolution with fc ov = const up to z max is 
assumed. 

• The sources arc assumed to be distributed homo- 
geneously in the universe, so that the propagation 
mode in extragalactic magnetic fields, from rectilin- 
ear to diffusive, does not influence the energy spec- 
trum at observation. The propagation is studied 



using a Monte Carlo approach with the inclusion 
of pair-production energy losses and nuclei photo- 
disintegration on CMB and EBL. 

At a first glance one may expect that the large num- 
ber of free parameters, such as generation index, parame- 
ters of cosmological evolution and coefficients of source nu- 
clei composition, can provide a broad variety of observed 
mass compositions and spectra. However, as it was demon- 
strated in pjj ]. the predictions are very much constrained 
due to photo-disintegration of nuclei on EBL and CMB 
radiations. The basic physics phenomena and their results 
are as follows. 

Gcncrically in the mixed models the mass composition 
becomes lighter at E > 10 EeV, because intermediate and 
heavy nuclei are destroyed by the EBL photons while pro- 
tons survive. In principle this situation may change only 
above 50 EeV, when GZK cutoff in the proton spectrum 
appears, while heavy nuclei, e.g. Iron, are still not photo- 
disintegrated by the CMB photons and may dominate. In 
realistic cases the dominant component in mixed models 
are protons. 

In the academic case of a fixed nuclei species at the 
source (Allard et al (2008) in j28|) secondary protons are 
the dominant component at all energies if the primaries are 
light or intermediate mass nuclei, and only in the case of 
Iron the protons are subdominant. It is interesting to note 
that because of an efficient destruction of light and inter- 
mediate nuclei the two-component model with only proton 
and Iron injected at sources gives a reasonable agreement 
with observations 



81 



In the mixed models the large primary proton com- 
ponent contributes more than the secondary protons, 
strengthening further effect of proton dominance. There- 
fore, the GZK feature is also present in the mixed- 
composition models. 

The proton dominance in mixed models is illustrated 
by Fig. [TTJ In the left panel decomposition of the diffuse 
flux in the elemental components is presented. The dom- 
inant component is primary protons, and the next sub- 
dominant is secondary protons and Helium. In the right 
panel the elongation curve (X max ) (E) is presented with 
data points from HiRes-MIA, HiRes-stereo and Fly's Eye 
superimposed. Starting from E > 3 EeV the elongation 
curve (X max ) (E) tends to a proton-dominated composi- 
tion in accordance with HiRes and Fly's Eye data. In 
contrast, the PAO data tend at the highest energies to an 
Iron-dominated composition. 

Transition from Galactic to extragalactic component in 
the mixed models depends on the choice of parameters. In 
most models transition occurs at the ankle, see Allard et 
al. in 28|. However, in the conceptuall y im portant paper 
by Allard, Olinto, Parizot (2007) from |28[ it was empha- 
sized that for strong source evolution and flat generation 
spectra the intersection of Galactic and extragalactic com- 
ponents occurs between 0.5 EeV and 1 EeV, i.e. at the 
second knee, as in the dip model. The transition in this 
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Figure 11. Left panel: Diffuse spectrum decomposed in its elemental components, calculated in the mixed-composition model by Allard, 
Olinto, Parizot (2007) from [28| . The proton dominance seen in the spectrum is a generic prediction of the mixed composition models. Right 
panel: Elongation curve (X max ) (E) in the same model as in the left panel. The mass composition evolves from almost a pure iron composition 
at E m 0.3 EeV to a lighter composition due to an enrichment by protons and light nuclei of extragalactic origin. At energy E a = 3 EeV 
the transition to pure extragalactic component is completed and chemical composition evolution proceeds further due to photo-disintegration 
of nuclei. At energy E ss 13 EeV all nuclei are disappearing faster than before and composition becomes strongly proton-dominated at 
E > 30 EeV. 



model is shown in Fig. [T2j The transition begins at 2 EeV 
and is accomplished at 3 EeV. 
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Figure 12. Diffuse propagated spectrum for mixed model by Allard 
ct al. (2007) [28ll with generation index -y g = 2.3 and cosmological 
evolution (1 + z) 3 up to 2 = 1.3 and with a frozen evolution at 
larger z up to z ma x = 6. The transition accomplishes at E Rj 3 EeV, 
i.e. below the ankle. At highest energies the spectrum is dominated 
by protons (see right panel of Fig. 111! and shows the GZK cutoff, 
compared in the figure with HiRes data. 

The dominance of protons (as an example see Fig. ITT1) 
was the reason why the (X max ) (E) predicted by mixed 
composition models is in a better agreement with HiRes 
data than with PAO data. The observations of PAO show 
that mass composition becomes heavier with increasing 
energy, and thus the existing calculations in the framework 
of mixed models agree better with the HiRes data. 

However, with the recent proposal (82^ . discussed in sub- 
section [331 the power of the mixed composition model for 
fitting the PAO data may change due to a possible en- 



hancement of the heavy nuclei production. 

5. Discussion and Conclusions 

The transition from Galactic to extragalactic CRs takes 
place in the energy interval (0.1 — 10) EeV. and it gives 
a key for understanding both galactic and extragalactic 
cosmic rays. 

In its low-energy part this range includes the high- 
energy tail of the Galactic CRs. The information obtained 
at these energies on the Galactic CRs involves the maxi- 
mum energy of acceleration, the chemical composition, the 
anisotropy and the mode of propagation in the Galaxy. It 
may unveil the whole picture of the origin of the Galactic 
CRs. 

The low-energy part of the extragalactic cosmic rays 
spectrum can give a key information on the existence of 
the pair-production dip and on the propagation of the CRs 
in extragalactic magnetic fields. 

Therefore, the experimental studies in the transition re- 
gion (0.1 — 10) EeV are of paramount importance in this 
field of research, with the mass composition measured by 
different methods being probably the most important task. 

There are four working detectors which cover par- 
tially the a bove -mentioned region: KASCADE-Grande 
IH, Tunka [l00|, Yakutsk ^ and IceTop/IceCube [lPl] . 
There are also projects to extend the observations of Tele- 
scope Array and PAO to low e nergi es, (0.1 — 1) EeV at 
TALE [lQ2j and at LE- Auger [l03j]. The Auger detec- 
tor has a great potential to explore the low-energy re- 
gion of the UHECR spectrum. At present there are al- 
ready two new detec tors at PAO collecting data at this 
energy; HEAT |l04^ . High Elevation Auger Telescope, 
detecting the fluorescent light at higher elevation an- 
gles; and AMIGA Auger Muons and Infill Ground 
Array, for the detection of the EAS muon component. 
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These detectors, together with TALE, Tvmka, Yakutsk 
and IceCube/IceTop will provide information on all ra- 
diations from EAS, including fluorescent and Cherenkov 
light, muons and radio radiation. One may expect that in 
this way the present controversy between the mass com- 
positions in the HiRes and Auger detectors will be unam- 
biguously solved. The recent measurement 76, 77J of the 
muon-production depth A^ ax and maximum zenith angle 
$max by on-ground Auger detectors is an important step 
in this direction. 

At the low-energy end of the UHECR the energy spec- 
tra are measured with an unprecedented accuracy for cos- 
mic ray physics. However, in fact even in experiments 
with the same technique, like HiRes and Auger, the en- 
ergy scales are different due to systematic errors (see left 
panel of Fig. [5]). Such difference is more natural for differ- 
ent techniques, e.g. scintillator and fluorescent ones in TA 
detector. The problem of the energy scale hopefully may 
be solved combining all observational methods. However, 
there is a physical 'standard candle' for the detector energy 
calibration, given by the fixed energy position of the pair- 
production dip. The recalibration factor k ca \ can be found 
by the spectrum shift to the energy at which the agreement 
between observed and predicted dips is the best [15j, [l_6 1 . If 
the observed dip is really caused by pair-production pro- 
cess, the fluxes measured in all experiments must coincide. 
This is what we see indeed in Fig. El in [HI M, [z3 and in 



Bliimer at al. 



Since the mass composition is expected to be different 
for Galactic and extragalactic CRs, its measurement gives 
a reliable way to test the transition. The best method at 
present to determine the CR chemical composition is given 
by measuring (A max ) (E). Unfortunately this method 
suffers from many uncertainties, including those in the 
value of fluorescent yield, absorption of UV light in the 
atmosphere and uncertainties in the interaction models, 
needed to convert the mass composition into the observed 
(A max ) (E). Systematic errors in the (X max ) measurement 
can be as large as 20 — 25 g/cm 2 , to be compared with 
the difference of about 100 g/cm 2 between the (A max ) (E) 
for protons and Iron. A better sensitivity to distinguish 
different nuclei is given by the width of the (A max ) (E) 
distribution, i.e. RMS(A max ) 

The most self-consistent conclusions on the nature of 
the observed UHECRs are obtained at present for HiRes 
detector: it has observed both signatures of proton propa- 
gating through CMB; namely pair-production dip and the 
GZK cutoff. Moreover, the conclusion about protons as 
primaries is confirmed by the direct measurements of pro- 
ton composition with the help of (A max ) (E) and RMS. 
GZK cutoff is also found in the integral spectrum with 
characteristic energy logE 1 / 2 = 19.73 ± 0.07, in excellent 
agreement with theoretical prediction logi^/2 = 19.72. 
Good agreement with the dip energy shape and with -E1/2 
allow to conclude that the HiRes basic energy scale is cor- 
rect. HiRes data agrees well with the dip model with AGN 
as the sources and with cosmological evolution of AGN 



taken from X-ray observations. 

The Auger data are quite different. Measurements of 
mass composition using (A max ) (E) and RMS and con- 
firmed by independent data of muon-production depth 
A^ ax and rise-time of the signal in water-Cerenkov tanks 
indicate nuclei as primaries with steadily increasing (A) 
starting from 4 EeV and up to 35 EeV, where mass com- 
position practically reaches pure Iron. The Auger mass 
composition, if confirmed, excludes the dip model and 
the observed high-energy steepening as GZK cutoff. At 
present there is no nuclei-based model which explains si- 
multaneously the Auger energy spectrum J{E) and elon- 
gation curve (A max ) (E). 

There are three models of transition: ankle, dip and 
mixed composition one. They differ most notably by the 
energy of transition and by the mass composition of the 
extragalactic component. 

The ankle model was the first one which appeared in lit- 
erature. It was proposed already in 1960s after the discov- 
ery of a flat feature, the ankle, at the Volcano Ranch array 
[20( 1 . The transition in this model is described in the most 
natural way, as the intersection of a steep (oc E~ 3A ) Galac- 
tic heavy-nuclei component (Iron) with a flat extragalac- 
tic proton component with a generation spectrum oc E~ 2 . 
The transition, as calculated in different models, occurs at 
E tr ~ (3 - 10) EeV. To fill the gap between 0.1 EeV (Iron 
knee) and E tv indicated above the ankle models need a new 
yet unknown high energy Galactic component beyond the 
Standard Model for Galactic Cosmic Rays. Recently the 
ankle model has been strongly disfavored by (A max ) mea- 
sured below the ankle (see right panel of the Fig. [TU)l due 
to small A max predicted by the model. It is also disfa- 
vored by the large (not observed) anisotropy predicted at 
E ~ (1 - 3) EeV pi. 

The dip model is based on the pair-production dip in- 
duced in the extragalactic spectrum of protons interacting 
with the CMB radiation. The dip is well confirmed by 
the observations of HiRes, TA, Yakutsk and AGASA (see 
Fig. |4]) . The successful energy calibration of all observa- 
tions (see Fig. [5] as an example) using the fixed energy 
position of the pair-production dip can be considered as 
another (indirect) evidence in favor of the dip model. 

The transition in the dip model occurs as an intersection 
of the steep Iron Galactic spectrum with a very flat extra- 
galactic spectrum at energy (0.1 — 1) EeV. The flatness 
of this spectrum can be provided as a ' mag netic flatness' 
in the case of diffusive propagation I90l |91| or due to the 
sources distribution over E max (83l . Il6| . A robust predic- 
tion of the dip model is the proton-dominated mass com- 
position above E > 1 EeV. It is confirmed by the HiRes 
and TA data; the calculated elongation curve (A max ) (E) 
is in good agreement with the bulk of data (see left panel 
of Fig. [TU)l . However, it strongly contradicts the highest 
energy Auger data, especially the RMS(A max ), presented 
in the right panel of Fig. [3 The Auger - HiRes controversy 
on the mass composition is a crucial point for the survival 
of the dip model. 
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In the mixed composition models the presence of various 
nuclei is assumed in the generation spectrum of the extra- 
galactic component. This assumption is very reasonable, 
because various nuclei are present in the gas where acceler- 
ation process operates. Mixed composition models with a 
wide choice of parameters are considered in [2j| ■ They in- 
clude different generation indices, cosmological evolution 
and most notably different relative abundances of nuclei 
at generation. However, as was explained in subsection 
14.31 a basic feature common to most mixed composition 
models is the proton dominance at the highest energies 
E > 10 EcV. This feature is in conflict with the Auger 
data at these energies that shows the mass composition 
getting progressively heavier with energy increasing and 
reaching Iron at E w 35 EeV. On the other hand, the 
mixed composition models fit reasonably well the HiRes 
and Fly's Eye data. The Galactic spectrum at the highest 
energies is assumed to be Iron dominated. The transition 
energy is a model-dependent quantity, being determined 
by 7 S and the parameters chosen for sources cosmolog- 
ical evolution. However, for the flat generation indices 
7 g = 2.1 — 2.3 and a reasonable choice of the cosmologi- 
cal evolution, the transition energy is about 3 EeV (as an 
example see Fig. [T2"]) . The elongation curve shown in the 
right panel of Fig. [TT] changes the slope at approximately 
the same energy, and the mass composition evolves to- 
wards the pure proton composition. The calculated elon- 
gation curve in Fig. [TT] agrees well above the transition 
with the data of HiRes stereo and Fly's Eye, though nowa- 
days the data of HiRes agrees better at E > 3 EeV with 
a pure proton composition. 

Much better quality of data on the (X max ) (E) is needed 
to distinguish the dip and mixed-composition models by 
I mM measurements. 

With the recent proposal [§2| for an enhancement 
of heavy nuclei production at the sources, the mixed- 
composition models may hopefully describe better the 
Auger elongation curve (X max ) (E) and its RMS. 

As short conclusion one may state that while the an- 
kle models of transition are severely disfavored, mostly 
by the (A max ) (E) data at E ~ 10 18 eV, the dip and 
mixed composition models are in a good agreement with 
the (A max ) (E) data from the Fly's Eye, HiRes and TA 
detectors. The present problem is the conflict of data on 
mass composition between the HiRes and TA on one side, 
and the PAO on the other. We do expect that problems 
with mass composition in the energy range (0.1—3) EeV, of 
fundamental importance to asses the transition issue, will 
be soon resolved with the HEAT, AMIGA, KASCADE-G, 
TALE, Tunka, Yakutsk and IceCube/IceTop experiments. 
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